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ABSTRACT 



The regularization of a new problem, namely the three-body problem, using 
'similar' coordinate system is proposed. For this purpose we use the relation 
of 'similarity', which ha s been introduced as an equivalence relation in a previ- 



ous paper (see 



Roman 



(1201 lh ). First we write the Hamiltonian function, the 
equations of motion in canonical form, and then using a generating function, 
we obtain the transformed equations of motion. After the coordinates transfor- 
mations, we introduce the fictitious time, to regularize the equations of motion. 
Explicit formulas are given for the regularization in the coordinate systems cen- 
tered in the more massive and the less massive star of the binary system. The 
'similar' polar angle's definition is introduced, in order to analyze the regulariza- 
tion's geometrical transformation. The effect of Levi-Civita's transformation is 
described in a geometrical manner. Using the resulted regularized equations, we 
analyze and compare these canonical equations numerically, for the Earth-Moon 
binary system. 



Subject headings: Restricted problems: restricted problem of three bodies . Stellar 
systems: binary stars. Methods: regularization 
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Introduction 



In a previous article (see iRoman I ( 1201 ll )). by introducing the "similarity" relation 
and applying it to the restricted three-body problem, the "similar" equations of motion 
were obtained. These equations were connected with the classical equation s of m otion 



Roman 



(1201 lh ). In 



by some coordinate transformation relations (see equations (17) in 
this paper were also defined 'similar' parameters and physical quantities, and 'similar' 
initial conditions and some trajectories of the test particles into the physical (x, Si,y) and 
respectively (x,S2,y) planes were ploted. 

Denoting Si and 5*2 the components of the binary system (whose masses are mi and 
777.2), the equations of motion of the test particle (in the frame of the restricted three-bo dy 



prob 



em) in the coordinate system (x, Si, y, z) are (see equations (11)-(13) in (IRoman 
201 lh ): 

d 2 x dii q x q(x — 1) 

X-- 1 



d 2 x 



2— 

dt 



1 + q (1 + q)r\ (1 + q)rl 



d 2 y dx 
— + 2 — 
dt 2 dt 



y - 



qy 



1 + q)r\ (1 + q)r% 



(2) 



where 



ri = a/ x 2 + y 2 + z 2 



q z 



d 2 z 



dt 2 (1 + q)r\ (1 + q)r. 



3 • 



r 2 



\/{x - l) 2 + y 2 + z 2 , q 



m 2 
mi 



(3) 



(4) 



In the 'similar' coordinate system (x', 5*2, y ' : z') t he equations of motion of the test 



particle are (see equations (14)-(16) in 




q'(x' - 1) 



d 2 x' ^dy' , 

Hi 2 + ~dT~ X ~l + q' (1 + q')r'? (1 + q')r' 2 3 



(5) 
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d 2 y' dx' 
dt 2 dt 



q' y' 



(1 + q')r? (1 + q')r'i 



(6) 



where 



d 2 z' 



q' z' 



dt 2 (1 + q'y? (1 + q')r' 2 



'x' 2 + y' 2 + z' 2 , r' = JJx' - l) 2 + y' 2 + ^ 2 , g' = — 

m 2 



It can be easily verify that the equations of coordinate transformation are: 



(7) 



x' — 1 — x , y' — y , z! = z . 



(9) 



One can observes that equations (l)-(3) and (5)-(7) have singularities in r x = 0, r 2 = 0, 
r[ = and r' 2 = 0. These situations correspond to collision of the test particle whith 5*i or 
5*2 in a straight line. The collision is due to the nature of the Newtonian gravitational force 
(F oc ^j). If the test particle approaches to one of the primaries very closely (r — > 0), then 
such an event produces large gravitational force (F — > 00) and sharp bends of the orbit. 
The removing of these singularities can be done by regularization. (Remark: the purpose of 
regularization is to obtain regular equations of motion, no regular solutions.) 



Euler seems to be the first (in 1765) to propose regular izing trans 



studying the motion of three bodies on a straight line (see 



Szcbchcly 



regularization metho d 



Celletti et al. 



(2011); 



ras become popul ar in recent years (see 



ormations when 



(1975)). The 



Jimenez-Perez et al. 



Waldvogel 



(2011); 



(120061 )) for long term studies of the motion of celestial 



bodies. These problems have a special merit in astronomy, because with their help we 
can studied more efficient the equations of motions with singularities. At the collision the 
equations of motion possess singularities. The problem of sing ularities plays an important 



role under computational, physical and conceptual aspects (see iMioc et al. 



(120021); 



Csillik 
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( 120031 )). The singularities occurring at collisions can be eliminated by the proper choice of 
the independent variable. The basic idea of regular izat ion procedure is to compensate for 
the infinite increase of the velocity at collision. For this reason, a new independent variable, 
fictitious time, is adopted. The corresponding equations of motion are regularized by two 
transformations: the time transformation and the coordinate transformation. The most 
important part of the regularization is the time transformation, when a new fictitious time 
is used, in order to slow the motion near the singularities. 



2. The 'similar' canonical equations of motion 



The regularization can be local or global. If a local regularization is done, then the 
time and the coordinates transformations eliminate only one of the two singularit i es. An 



example for the local regularization is the Birkhoff's transformatio n (see 



The global re gularization eliminates both singularities at once (see 



Birkhoff 



Castilho et al. 



915b). 



Csillik 



(119991 ): 



( 120031 )). Because our singularities are given in terms of — , — , — , — , in this paper a 



global regularization will be done. 

In order to do this, we need to replace the cartesian equations (l)-(3) and (5)-(7) with 
the corresponding canonical equations of motion. The canonical coordinates are formed 
by generalized coordinates 91,92,93 and generalized momenta Pi,P2,P3- The Hamiltonian, 
defined by the equation: 

3 „ 3 



dC 



i=l 



(10) 



1=1 



Boccaletti et al. 



will becomes (see 
coordinates system rotates): 

1 



(119961 1 p. 266, for the generalized momenta, when the 



M = ^(P? +P2 +P3) +Pi^2 - 9iP2 + y + y -^(91,92,93) • 



Til 
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Here 



with 



^(91,92,93) 



9i 



l + q 



4 



+ 



2q 



1 + q)r 1 (1 + q)r 2 



n \/ q\ + ql + qh r 2 = \l (91 - l ) 2 + 4 + 4 ■ 

Here the generalized coordinates and the generalized momenta were: 

qi = x , q 2 = y , 93 = z , Pi = q\ - q 2 P2 = 92 + 9i , Ps = 4a- 



Then, the canonical equations 



Qi 



m 

dpi 



Pi 



on 

dqi 



€{1,2,3} 



have, in the (91, «Si, 92, 93) coordinate system, the explicit forme: 

Pi + 92 
P2 ~ 9i 

P3 
P2 ~ 

~Pl 



dqi 

dt 
dq 2 

dt 
dq 3 

dt 
dpi 

dt 
dp 2 

dt 
dp 3 

dt 



9 

l + q 
1 



1 . Si _ 9 

l + q r\ l + q 

92 9 92 



9i 



l + q rf l + q r| 



93 



9 



93 



(12) 
(13) 

(14) 

(15) 

(16) 
(17) 
(18) 
(19) 

(20) 

(21) 



l + q r{ l + q r\ 

It is easy to verify that using the relations (14), the explicit canonical equations become the 
cartesian equations (l)-(3). 

In order to write the canonical equations in the 'similar' coordinate system 



[gi^S 2 ,q 2s ,q3 S ), we have in view the theoretical considerations from the article (IRoman 



sec 



2011). The index s refers to 'similar' quantities. Then, the 'similar' Hamiltonian will be: 



Boccaletti et al. 



( 11996( 1 p. 266): 



Us = „ (Pi, + P 2s + P 3s ) ~ (Pls92s - qisP2s) + -7T + ~7T ~ ^e(,QU, 92s, 93 



(22) 
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where 



qu 



1 + q' 



2 ^ 2 2g' 

+ + 



1 + q')r is (1 + q')r 2s 



with 



n, = \ qls + qls + qL r 2s = \ (qu - 1) 2 + qls + ql 



Here the generalized coordinates and the generalized momenta were: 



q ls = 1 - q x , q 2s = <?2 , q-3s = <?3 , Pis = ~Pl P2s = P2 ~ 1 , P3s = P3- 



Then, the canonical equations 



q%. 



Pis 



dm 



, 1 €{1,2,3} 



r\ 1 fis r% 

opis oqi 

have, in the (qu, S 2 , q 2s , <?3 S ) coordinate system, the explicit forme: 
dqi s 



dt 
dq 2 s 

dt 
dq3s 

dt 
dpu 

dt 

dp2s 
dt 

dp3s 
dt 



Pu ~ <?2 

P2s + qi 

P3s 

-P2s - 
Pis 



qis 



q qu 



l + q' l + q' r\ s 1 + q' 

1 q2s q' q2s 



2s 



1 + q' rf s 1 + q' r% s 



q3s 



q qss 



(23) 



(24) 



(25) 



(26) 



(27) 
(28) 
(29) 
(30) 

(31) 

(32) 



1 + q' rf s 1 + q' r| s 

It is easy to verify that using the relations (27), (28), (29), the explicit canonical equations 
(30), (31), (32), become the cartesian equations (5)-(7). 

Remark: From e quations (13) and (24) it is easy to observe that r\ = r 2s and r 2 = ri. s 



(see also Figure 1 in (IRoman 



201l|)). 
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3. Coordinate transformation 

The equations of motion (19)-(21) and (30)-(32) have singularities in n and r 2 , 
respectively in r ls and r 2s . We shall remove these singularities by regularization. Several 
regularizing methods are known (see Stiefel et al. 1971). In this paper we shall use the 
Levi-Civita's method, applied when the bodies are moving on a plane. The two steps 
performed in the process of regularization of the restricted problem are the introduction 
of new coordinates and the transformation of time. The combination of the coordinate 
(dependent variable) transformation and the time (independent variable) transformation 
have an analytical importance and increase the numerical accuracy. For simplicity we shall 
consider that the third body moves into the orbital plane. 



3.1. Case 1 - coordinate transformation in the coordinate system with origin 

in Si 



For the regularization of the equations of motion in the (qi, Si, q 2 ) coordinate system, 
we shall introduce new vari ables Qi and co nected with the coordinates qi and q 2 by the 



relations of Levi-Civita (see 



Levi-Civita 



(119061 )): 



qi = Ql-Ql, Q2 = 2QiQ 2 



Let introduce the generating function S (see 



Stiefel et al. 



jl97lh . p.196): 



S = -pif(Qi,Q2) -P29(Qi,Q2) , 

a twice continuously different iable function. Here / and g are harmonic conjugated 
functions, with the property 

df dg 



(33) 



(34) 



dQi 
df_ 
dQ 2 



dQ 2 

dg_ 

dQi 
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The generating equations are 



Qi = — 



dS 



Pi = 



dS 



' dQi ' 

with Pi, P2 as new generalized momenta, or explicitly 

dS 



IE {1,2}, 



(35) 



Qi 
Q2 
Pi 
P2 



dpi 

ds 

dp 2 
dS 

dS 

W2 



df dg 



Pi 



dQ 1 
df 



+ P2 



dQi 
dg 



where 



an 



dQ 2 " dQ 2 
df dg 



Pian + p 2 a l2 



(36) 



Ol2 = 

Let introduce the following notation: 

an ai2 



dQ l dQ 2 
df dg 



— 0-12 On 



Pi 



,P2 



P 2 



dQ 2 dQ 1 



, D = detA = a n + a 2 2 , 



A T 



, p = ^-p; p = — P,pi+ P 2 2 = (P 1 2 + P 2 2 )/D, 



(37) 



where A T represents the transpose of matrix A. The new Hamiltonian with the generalized 
coordinates Q\ and Q 2 and generalized momenta Pi and P2 is: 



H(Q 1 ,Q 2 ,Pi,P 2 ) 



1 
2D 
1 



+ 



l + q 



f~ 



1 



q 



l + q r x l + q r 2 2(1 + q)' 



(38) 



where r x = \J f 2 + g 2 , r 2 = \/(J~- l) 2 + g 2 and D = A{Q\ + Q 2 ) and the explicit canonical 
equations of motion in new variables become: 



dQi 
dt 



1 
2D 



2Pl + 4 (/2+92) 
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dQ 2 
dt 
dP l 
~df 



dt 



1 
2D 
Pi 



2Po - 



d_ 

dQ 



(f + 9 2 ) 







2D dQ 1 dQ 2 



+ 







l + q dQ 1 
Pi d 



2D dQ 2 dQ 2 



+ 



d 



1 



l + q dQ 2 \r 1 

Using Levi-Civita's transformation / 
(33)), the equations (39) becomes: 



+ 



P2 


d 


2D ' 


dQxdQx 


d 




dQi 


(Q 


P2 


d 


2D ' 


dQ 2 dQ l 


d 




dQ 2 


(k) 



(f 2 + 9 2 )- 



g df 
l + qdQ 1 



if + 9 2 ) 



l + qdQ 2 



+ 



(39) 



Qi — Qi - Q\ > 9 — 92 = ^QiQi (see relations 



dQ l 


-h 


Q2 


dt 




2 


dQ 2 


P2 


Qi 


dt 


D 


2 


dP l 


P2 


2qQi 


dt 


2 


l + q 


dP 2 


Pi 


2qQ2 

l + q 


dt 


2 



2 Q, 2q Q 1 (f 1 -1) , (Pf + P^Q, 



■ q r\ l + q 
2 Q 2 



^3 



+ 



4f 2 



(40) 



2q g 2 (r 1 + i) + (p 1 2 + p 2 2 )g 2 



l + q 



^3 



4r 2 



where r 1 =Q\ + Q 2 , r 2 = y/(Q\ -Q\- l) 2 + 4Q 2 Q 2 , 
with the new Hamiltonian 

Ksi = fi + + ^(PiQ 2 - P2Q1) + Q 



8(Q? + Qi) 
1 1 



1 + g 



1 + 9 Q 2 + Ql l + q ^(Ql-Ql-lf + AQ\Ql 2{l + qf 



(41) 



3.2. Case 2 - coordinate transformation in the 'similar' coordinate system 

For the coordinate transformation in the (qi s , S2, q2s) coordinate system, we introduce 
the generating function S s in the plane (q sl , S 2 , q S 2), hi the following form 



<S S — —Pslfs(Qsl,Qs2) — Ps29s(Qsl, Qs2) , 



(42) 
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where f s and g s are harmonic conjugated functions. The generating equations are 



or explicitly 



q s i 

Psl 
Ps2 



P* 



dS s 

dpsi 
dS s 

dS s 
dQsi 



dpsi 
dSs 
dQ s 



ie{l,2} 



fs(Qsl, Qs2) 
9s(Qsl, Qs2) 

df s . dg 

= Psl^—+Ps2 



dQ 



dQsi dQsi 
dfs . dg s 

PslT^ \~Ps2- 



s2 



OQ 



»2 



'dQ 



Psihi + p s2 bi2 

~Pslbl2 + Ps2bll 



s2 



where 



611 
612 



dfs 


dg s 


dQsi 


dQ S 2 


dfs 


dg s 


90 »? 





Let introduce the following notation, (jSzebehely 

611 612 



(J1967I), p. 373) 



B = 



,D S = detB = b\ x + b\ 2 , 



H2 o u 



Psl 
Ps2 



The new Hamiltonian for the case 2 may be written 
1 



V Ps 2 = ^Ps 2 , P 2 sl + P 2 s2 = (Psl + Ps2)/Ds 



H 



■S2 



IDs 
1 



Psl + P S 2 ~ Psl 



J^ifs 2 + 9 2 s ) + Ps2^(f 2 s + ft 



+ 



1 + q 



,1 f s 



1 



q' 



q 



/2 



(43) 



(44) 



(45) 



_ (46) 

1 + q' r.i 1 + q' f s2 2(1 + q'f V ; 

where r si = y/ f% + g%, r S 2 = a/ (f s — l) 2 + g 2 s and D s = 4(Q 2 1 + Q 2 s2 ) and the Hamiltonian 

equations in new variables become 



dQ 
dt 



si 



2D, 



2P s i- 



_d 

dQ 



s2 



ifs+gl 
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dQ s2 


1 


dt 


2D S 


dP s i 


Pel 


dt 


2D S 




1 



2P. 



s2 



g 



si 



d 



+ 



!s2 



s2 



d 



2D S 8Q sl dQ sl 



Ul + al)- 



d 



dP. 



s2 



l + g' dQ sl 

Pel 9 



Q 



d 



1 + q' dQ sl 

Ps2 



1 

d 



dt 



2D,. 0Q, 2 0Q s2 (/;> + ^ 2D S ' dQ s2 dQ sl (/ * + & 



q' df a 
1 + q' dQ sl 



q' df, 
1 + q' dQ s2 



+ 



+ 



+ 



d 



q' 



d 



(47) 



1 + q' dQ s2 \r sl J 1 + q' dQ s2 
Because the singularity of the problem is given by the terms ljr s \ and l/r s2 , we will made 
a global regularization using the Levi-Civita's transformation 



fs — Qsl Qs2i 

The 'similar' Hamiltonian equations are given by 



9s — 2QslQs2 



(48) 



dQsi 


P.i 


Qs2 


dt 


D s 


2 


dQs2 


P S 2 
= — + 

D s 


Qsi 


dt 


2 


dP sl 


Ps2 


Iq'Qsi 


dt 


2 


1 + q' 


dP s2 




2q'Qs2 


dt 




1 + q' 



2 Q 



si 



1 + q' r 2 



2q' Q s i(f sl -1) + {n + Pl)Qsi 



si 



1 + q' 



r 3 s2 



4rj 



2 Qs2 2q' Q s2 (f sl + 1) , (P? 1 + P s 2 2 )Q s2 



1 + q' r 2 



si 



—3 

r S 2 



+ 



4rf 



(49) 



where r sl = Q 2 sl + Q 2 s2 , f s2 = ^(Q'j, - Q'j 2 - l) 2 + 4Q 2 sl Q 2 s2 , 



with the new Hamiltonian 

p2 _|_ p2 j q/ 

U * 2 = S (Q 2 sl + Q 2 s2 ) + 2 {Ps2Qsl ~ PslQs2) + TT?^ 1 " ^ " 



1 



1 



q' 



1 



J2 



i + q' Q 2 si + Q 2 S 2 i + q' V(Q 2 si-Q 2 s2-i) 2 + iQ 2 iQ 2 s2 2 (! + ^ 



(50) 



4. Time transformation 



The transformation of the independent variable is necessary to achieve regularization. 
It is a slow-down treatment of the physical problem, a new time scale in which the motion 



13 



slows down (IMikkola et al. 1119961 ) 



4.1. Case 1 - time transformation in the coordinate system with origin in Si 



To resolve t 



le Hamiltonian equations fl4"0l). we introduce the fictitious time r 



see 



Szebehely 



019671); 



Waldvogel 



(11972 



19821 ): 



Erdi 



( 120041 )). and making the time 



transformation ^ = rfrf, the new regular equations of motion are 



dQi 
dr 


Pinrl 

4 


Q 2 r 2 r\ 
2 




dQ 2 


P 2 rir\ 


Qir\rl 




dr 


4 


2 




dPi 
dr 


P 2 r\r\ 
2 


2 ^! f 2-3 

l + q 1 2 






2Qif\ 


2gQ 1 (r 1 - l)r? 


{Pi + PDQirl 




1 + q 


1 + g 


4 



dr 



Pif{ri , 2gQ 2 _ 2 



+ 



r r 



3 



1' 2 



2Q 2 r 2 i 2qQ 2 (fi + l)r\ (P 2 + P 2 )Q 2 r\ 



l+q l+q 
The explicit equations of motion may be written 

Pi 



d 2 Qi 



dr 2 



d 2 Q 2 
dr 2 



ldPi_ _ 3 lrfQ 2 _ 2 _3 . 

; — T\T H ; — ^1^9 + . 

A dr 2 2 dr 1 2 \ 2 



dQ 



+ 3 | ^ + Q 2 r^ f (Q? + Q1Q2 " Qx, dr 



ldP 2 o ldQi_ 2 _ 3 

rir 2 — — — — 7*1 r 



+ 3 



4 dr 
P2 



v 2 



2 



2 rfr 

f (Ql + QiQl - Qi 



(Q2 + QIQ2 + Q2 

dQi dQ 2 
VI— ; r V2 



dr 



dr 



dr 



+ (Ql + QIQ2 + Q2 



r 2 + 

dQ 2 
dr 



r\ + 



dQ 2 
' dr 



(51) 



rir 2 



rir 2 (52) 



Remark: It is easy to see that now, the equations of motion have no singularities. 
For the application of the above problem in a binary system, we can obtain the solution 
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in the form 



qi (t) = Qi(t)-Q 2 2 (t) 
q 2 (t) = 2Q 1 (t)Q 2 (t) 



(53) 



4.2. Case 2 - time transformation in the 'similar' coordinate system 

Introducing the fictitious time r and making the time transformation ^ = r^f^ 2 , the 
new regular equations of motion are obtained in the form 

dQsi Psir s ir 3 s2 _ Qs2r 2 sl f 3 s2 

dr 4 2 

dQs2 Ps2r s ir 3 s2 Qsir 2 sl f 3 s2 

dr A 2 

dP s i Ps2f 2 sl f 3 s2 _ 2q'Q sl _ 2 _ 3 

dr 2 1 + q' TslTs2 

2Q sl f 3 s2 2q'Q sl (f sl -l)f 2 sl , (P^ + P^Qsifl, 



1 + q' 1 + q' 



+ 



dPs2 _ Psir 2 sl rj 2 2q'Q s2 _ 2 _ 3 
dr ~ + ' 2 + l + q'^ 1 ^ 2 

2Q s2 fl 2 2 q 'Q s2 (r sl + iy sl {P 2 1 + P 2 2 )Q S 2f 3 s2 
1 + q' 1 + q' 4 ' 1 ' 

The explicit equations of motion are given by 

d 2 Q s \ _ ldP s \_ _ 3 1 dQ S 2_ 2 _ 3 / P s i 0/n _ \ c?Qsi „ dQ s2 \ _ 3 



+3 (^f - Q, 2 r,ij ((<& + e,i<3= 2 - Q.i)^? 1 + (<& + QliQ,2 + 

d 2 Qs2 ldP s 2_ _3 lc?Qsl_2 -3 / Ps2 _ \ / „ G?Qsl „ 



<^3s2 \ _ _ 



dr2 4 rfr ^ 2 + + {-f + 2 Qsl r sl j \q a ^ + Qs2 ^ j n 2 + 

'Ps2 , ^ _ \ /^ 3 , ^ 



+3 I -f + Q sl f sl \ UQ 3 sl + Q sl Qj 2 - Qsi)-^ + (Q 6 s2 + Q 2 sl Q S 2 + Qs2)~^ ) f sl r s2 (55) 



Remark: It is easy to see that the 'similar' equations of motion have no singularities. 
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For the application of the above problem in a binary system, we can obtain the solution 
in the form 



Qsi(t) = Q 2 sl (t) - Q 2 s2 (t) 
q s2 (t) = 2Q sl (t)Q s2 (t) 



(56) 



5. Numerical experiments 



For the numerical integration (Ea rth-Moon bina ry system), considering that the third 



body moves into the orbital plane (see 



Kopal 



(119781 )). we used the initial values: 



gio = 0.6, q 20 = 0.4, p w = 0.1, p 20 = 0.6, t £ [0, 2tt], q = 0.0123 . 

For the numerical integration (Earth- Moon binary system) in the "similar" coordinate 
system we use the initial values (see eqs. (25)): 

q 10s = 1.6, q 20s = 0.4, p 10s = -0.1, p 20s = -1.6, te[0,2vr], q' = 81.30 . 

For the numerical integration (Earth-Moon binary system) in the regularized coordinate 
system (equations (52)), we use the initial values (see also eq. (33)): 

Qio = 0.813, Q 20 = 0.246, P 10 = 0.458, P 20 = 0.926, r€[0,2vr], q = 0.0123 , 

and in the 'similar' regularized coordinate system (equations (54)): 

Q 10s = 1.275, Q 20s = 0.157, P 10s = -0.757, P 20s = -4.047, r e [0, 2tt], q' = 81.30 . 



5.1. Considerations on the initial conditions 



In Figure 1 we can compare the trajectories of the test particle in the coordinate 
systems with origin in Si (figures a, c, e), and S2 (figures b, d, f). The point Pi correspond 
to the initial conditions. 
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We consider the trajectories given in (IRoman Il201ll ) in figure 6 (in the coordinate 
systems (x,Sx,y) and (x',S2,y')) and we represented them in the coordinate systems 
(ql, Si, q2) and (qls, S 2 , q2s) (see Figure 1 a and b). In this purpose we obtained the initial 
conditions as follows: 



<?io = Xo = 0.6 ; q 20 = yo = 0.4 ; <ji = v 0x = 0.5 ; q 20 = v 0y = 



and from eqs. (16)-(17): 



P10 = <?io - <?2o = 0.1 ; P20 = <l2o + qio = 0.6 



and 



qios = x'o = 1-6 ; q20s = y' = 0.4 ; q 10s 
and from eqs. (27)- (28): 



J 0x 



-0.5 ; q 2 os = v' 0y = , 



Pws = qws + q20s = -0.1 ; p 2 o s = <?20s - qws = -1-6 . 

In order to obtain the initial conditions, when we make the coordinate transformation, we 
solve the systems: 



Qio — Qio ~ Q20 I -fio — 2pio<5io + 2p2oQ 



20 



Q20 — 2Q10Q20 ^ -P20 — — 2P10Q20 + 2p2oQio 

(see eqs. (33) and (36)) for the trajectory in (Qi, 5i, Q2) coordinate system (Figure lc) and 



QlOs — QlOs Q20S 



PlOs — 2pi0sQl0s + 2p20sQ20s 



Q20s — 2Ql0sQ20s ^ P20s — — 2pi0sQ20s + 2p20sQl0s 

(see eqs. (44) for the trajectory in (Qi s , S 2 , Q2s) coordinate system (Figure Id). 

Obviously, the initial conditions remain the same if we change the real time t to the 
fictitious time r, but the motion is slowed. In Figure le and Figure If we represented the 
motion in real time t with thin line and the slowed motion with thick line (corresponding 
to the same period of time). 
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5.2. Considerations on the geometrical transformation 

Let us analyze the Figures la and lc. For this purpose we consider a point A(qi,q 2 ) 
on the graphic show in Figure la, and B(Qi,Q 2 ) its corresponding point in Figure lc. We 
have (see Figure 2a and 2b): 

+ ? 2 2 QiQz 2tanBS 1 Q 1 -^y- 
taniqxbxA) = — = —k — , = = tan(2 BSiQi) 

Qi Ql-Ql l-tan 2 BS 1 Q 1 

and it results: ASiqi = 2 BS\Q\. We used the counterclockwise directions for measuring 
the angles. 

The Levi-Civita geomet rical transformation originate in the conformal transformation 



see 



Boccaletti et al. 



( 11996( 1. p.164): 

z = q x + i q 2 = (Qi + i Q2) 2 

where (q±, Si, q 2 ) is the physical plane and (Qi, Si, Q 2 ) is the parametric plane. From this 
relation we have: qi — Q\ — Q 2 , q-i — 2 QiQ 2 , and \SiA\ = \SiB\ 2 . It means that the 
geometrical transformation squares the distances from the origin and doubled the polar 
angles. 

If, having the trajectory in the physical plane, we want to draw the trajectory into 
the parametric plane, we have to choose a point Ai on the trajectory in (qi,Si,q 2 ) plane, 
measure the angle AiSiqi and the distance SiA^, and then draw a half-line B[Si in the 
(Q1S1Q2) plane, so as A^Siqi = 2 B[SiQi. O n this half-line, we have to measure the 
distance SiB{ = y/ (S1A4), and obtain the point Sj. Than we have to repeat the procedure 



for % = l,n, n € N. Of course the computer will do this better and faster than we can do 
it, but the above considerations help us to understand what it happened. 

The vertex of the polar angles have to be centered into the more massive star, so the 
angles qiAq 2 and qi s A s q 2s and respectively QiBQ 2 and Qi s B s Q 2s are 'similar' polar angles. 
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So, if we intend to study the regularization of the circular restricted three-body problem 
using 'similar' coordinate systems, we have to add to 'similar' parameters postulated in 
section 3 in (Roman, 2011), the 'similar 7 polar angles, measured between the abscissa and 
the half-line passing through the center of the most massive star and the test particle. 

5.3. The effect of the Levi-Civita's regularization 

In order to see what is the effect of geometrical transformation, let us analyze the 
graphics from Figure 3. In Figure 3a there are represented some circles; their equations 
are: q\ + q\ = r 2 , where r G {0.2; 0.4; 0.6; 0.8; 1; 1.2; 1.4; 1.6; 1.8; 2}. In figure 3b there are 
represented the circles having the equations Q\ + Q\ = u 2 , where u = y/r, like geometrical 
transformation of Levi-Civita's regularization postulated. One can see that the circles in 
Figure 3b go away from the center and draw near the circle having radius u — 1. If in the 
center of the circles there is a problem (a singularity), it can be easily examined. 

In Figure 3c and Figure 3d there are represented some half-lines, having equations: 
q 2 = m qi, respectively Q 2 = n Q\, where m G {0; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8}, and 
n = 2 m, like geometrical transformation of Levi-Civita's regularization postulated. One 
can see that the half-lines in Figure 3d go away from the abscissa's axis. If there is a 
problem (a singularity) on the abscissa's axis, it can be easily examined. 

There are only two points invariant with respect to the geometrical transformation of 
Levi-Civita's regularization: Si(0;0), and S^ljO), respectively in the 'similar' coordinate 
systems Si(l; 0), and S 2 (0; 0). Then, the geometrical transformation go away the trajectory 
from the points where there are singularities. 

In what concern the time transformation, as one can see in Figure le and If, the role 
of this transformation is to slow-down the motion of the test particle. With thin line is 




Fig. 1. — Trajectories in different coordinate systems with origin in SI and S2 
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represented the trajectory of the test particle when the time integration is 2 tt, and with 
thick line is represented the trajectory of the test particle when the time integration is 40 tt. 
As one can see, after 40 tt we are still far away from the point where it is possible to have a 
singularity, if the coordinate system has the origin in Si, but not so far away if the origin of 
the coordinate system is in S2. 



6. Concluding remarks 



This paper continue the study of the relation of 'similarity', postulated in (IRoman 



201ll ) . by applying it to the Levi-Civita's regularization of the motion's equations of the 
test particle, in the circular, restricted three-body problem. Many papers in the last decade 
have studied the restricted three-body system in a phase space. During these studies, 
difficulties have arisen when the system approaches a close encounter. 

Using the regularization method in the 'similar' coordinates system, we give explicitly 
equations of motion for the test particle. We study numerically the regular equations of 
motion, we written in canonical form, and obtained that the integrator using regularized 
equations of motion are more efficient. The 'similar' Hamiltonian (see eq. (22)) give us the 
'similar' canonical equations (27)-(32), which have some different signes than the canonical 
equation (16)-(21). The coordinate transformation used in the Levi-Civita's regularization 
create a new form of the 'similar' Hamiltonian's equations (eqs. (49)). Finally, the time 
transformation used in the Levi-Civita's regularization gives us the regularized equations of 
motion (51) in the coordinate system with origin in Si, and (54) in the 'similar' coordinate 
system. 

In order to explain the shape of the trajectories in a concrete example, the 'similar' 
polar angle is introduced. 
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Our method may provide new directions for studies of circular restricted three-body 
integration using similar coordinate systems. It is an important tool for developing efficient 
numerical algorithms. 
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Fig. 2. — How to obtain the B point in (Q1,S1,Q2) plane from the A point from (ql,Sl,q2) 
plane 



- 23 - 




-24 - 

REFERENCES 

Birkhoff, G.D.: Rend. Cir. Mat. Palermo, 39,1 (1915) 

Boccaletti, & D., Pucacco, G.: Theory of orbits, vol.1, Springer- Verlag Berlin Heidelberg 
New York (1996) 

Castilho, C, & Vidal, C: Qual. Theory Dyn. Syst., 1, 1 (1999) 

Celletti, A., Stefanelli, L., Lega, E., & Froeschle, C: Celest. Mech. Dyn. Astron., 109, 265 
(2011) 

Csillik, L: Regularization methods in celestial mechanics, Cluj: House of the Book of 
Science (2003) 

Erdi, B.: Celest. Mech. Dyn. Astron., 90, 35 (2004) 

Jimenez-Perez, H., & Lacomba, E.: J. Phys. A, 44, 265 (2011) 

Kopal, Z.: Dynamics of Close Binary Systems, Reidel, Dordrecht (1978) 

Levi-Civita, T.: Acta Math. 30, 305-327 (1906) 

Mikkola, S., k Aarseth, S.J.: Celest. Mech. Dyn. Astron., 64, 197 (1996) 

Mioc, V., & Csillik, L: RoAJ, 12, 167 (2002) 

Roman, R.: Ap&SS, DOI: 10.1007/sl0509-011-0747-l (2011) 

Stiefel, L., & Scheifele, C: Linear and regular celestial mechanics, Berlin: Springer (1971) 
Szebehely, V.: Theory of orbits. Academic Press, New York (1967) 

Szebehely, V.: Regularization in celestial mechanics, in Lecture Notes in Mathematics, vol. 
461, p.257-263 (1975) 



Waldvogel, J.: Celest. Mech. 
Waldvogel, J.: Celest. Mech. 
Waldvogel, J.: Celest. Mech. 



-25 - 

Dyn. Astron., 6, 221 (1972) 
Dyn. Astron., 28, 69 (1982) 
Dyn. Astron., 95, 201 (2006) 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



